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We investigate the mixed state properties in a type II multiband superconductor with uniaxial 
anisotropy under the Pauli paramagnetic effects. Eilenberger theory extended to a multiband su¬ 
perconductor is utilized to describe the detailed vortex lattice properties, such as the flux line form 
factors, the vortex lattice anisotropy and magnetic torques. We apply this theory to SrjRuCU to 
analyze those physical quantities obtained experimentally, focusing on the interplay between the 
strong two-dimensional anisotropy and the Pauli paramagnetic effects. This study allows us to 
understand the origin of the disparity between the vortex lattice anisotropy (~60) and the H c 2 
anisotropy (~20). Among the three bands; 7 with the effective mass anisotropy ~180, a with ~ 120, 
and f3 with ~60, the last one is found to be the major band, responsible for various magnetic re¬ 
sponses while the minor 7 band plays an important role in the vortex formation. Namely, in a field 
orientation slightly tilted away from the two dimensional basal plane those two bands cooperatively 
form the optimal vortex anisotropy which exceeds that given by the effective mass formula with 
infinite anisotropy. This is observed by small angle neutron scattering experiments on Sr 2 Ru 04 . 
The pairing symmetry of S^RuCU realized is either spin singlet or spin triplet with the d-vector 
strongly locked in the basal plane. The gap structure is that the major /3 band has a full gap and 
the minor 7 band has a d x 2_ y 2 like gap. 

PACS numbers: 74.25.Uv, 74.70.Pq, 74.25.Ha, 61.05.fg 


I. INTRODUCTION 

It is now widely recognized that multiband supercon¬ 
ductors are omnipresent^. This recognition may be trig¬ 
gered by MgB^ where there exist distinctive two bands; 
the 3D 7 r-band and 2D like cr-band^. They play differ¬ 
ent role in forming superconductivity, in particular in 
magnetic properties under an applied field, such as sym¬ 
metry of vortex lattices^ or the form factors probed by 
small angle neutron scattering (SANS) experiments^. To 
understand its detailed magnetic response, a two band 
model is indispensable. In fact different dimensionality 
of the band structures between the 7r-band and tr-band 
gives rise to rotation of the triangular vortex lattice un¬ 
der varying field^. The form factors of SANS experi¬ 
ments clearly demonstrate a gradual change of the two 
components of the 7r-band and cr-band as field varies^. 
It is also true for other materials among unconventional 
and conventional superconductors where multiband de¬ 
scription is essential, such as heavy fermion superconduc¬ 
tors^ and iron pnictides^. 

We have been seeing that the Pauli paramagnetic ef¬ 
fect (PPE) is important when combined with this multi¬ 
band effect in certain superconductors, which give rise to 
a variety of unexpected phenomena. Typical examples 
are the oldest heavy fermion superconductors CeC^Si^ 
and UBei^i, and also KFe 2 As 2 & 2 , which necessitate the 
multiband description in fully understanding of their vor¬ 
tex properties. Those include the hidden first order tran¬ 
sition phenomenon^ and the disparity^ between the vor¬ 
tex lattice anisotropy and Fermi velocity anisotropy as 
discussed later. 


Here we study the interplay between PPE and multi¬ 
band effects in a uniaxial anisotropic superconductor in 
which each band has a different uniaxial anisotropy: For 
the case of two bands which we consider in this paper we 
can envisage two possible situations as schematically il¬ 
lustrated in Fig. [Q Let us consider the first case shown in 
Figs.[T|(a) and (b). In the absence of PPE the two orbital 
limited upper critical fields and H° cross at A in 

H versus ff plane (f l is the angle from the ab plane), each 
of which is characterized by the effective mass anisotropy 
r, for two bands i = 7 and /?, assuming F 7 > Tp. As 
indicated in Fig.[]](a) the four divided regions are char¬ 
acterized by each T;. In particular, along H C 2 (Ul) the 
characteristic anisotropy of the total system is switched 
at the intersecting point A from T 7 to as SI increases. 
When traversing at a higher H , the anisotropy T 7 for the 
7 band is sensed only. 

Now let us switch on PPE, then both orbital limited 
H° 2 h (Ul) are suppressed towards lower fields, especially 
r 7 if we assume that the superconducting gaps for the 
two bands such that Ap > A 7 . This is because the Pauli 
limited fields H? > H^. The resulting phase diagram is 
shown in Fig. [I] (b). As seen from it the crossing point 
A is removed and three regions are now occupied by F^ 
while the T 7 region is hidden deep inside at lower fields 
and finite fFs. In particular, along H C 2 (Ul) the region 
persists all the way from Q = 0° to fl = 90°. Thus 
the higher field scan is sensing only the T p anisotropy 
while a lower field scan is sensing F^ —>■ F 7 —F^ as fi 
increases. This non-trivial anisotropy evolution is caused 
by the interplay between the effective mass anisotropy T 
and PPE. 
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FIG. 1: (Color online) Schematic phase diagrams in the H 
versus Q, plane. The /3 band major scenario without PPE 
(a) and with PPE (b). The 7 band major scenario without 
PPE (c) and with PPE (d). The effective mass anisotropy is 
assumed to be T 7 > P^. 


Another possible phase diagrams are depicted in 
Figs. Etc) and (d) where < A 7 is assumed, keeping 
r 7 > r^. The 7 (/ 3 ) band with A 7 (A^) is now major 
(minor). In the absence of PPE shown in Fig. (He) the 
two orbital limited curves are not crossed, thus 

the higher (lower) field region is occupied by r 7 (r i a). In 
the presence of PPE those curves are both suppressed 
downwards. Thus the two regions (r 7 and T /3 ) are sim¬ 
ply shifted downwards, keeping its phase diagram topo¬ 
logically unchanged. Note that H& < H£. Therefore, 
a higher (lower) field scan is exclusively sensing the T 7 
(r^) anisotropy of the total system as f l increases. 

The purpose of this paper is to investigate the physics 
of Sr 2 Ru 04 through the studies of the mixed state prop¬ 
erties in this multiband superconductor where the in¬ 
terplay between the effective mass (or Fermi velocity) 
anisotropy and PPE is important in understanding of 
the vortex lattice state in F vs O plane. This su¬ 
perconductor is known to have the strong two dimen¬ 
sional uniaxial anisotropies for the three bands!!, the 
7 band (r 7 =180) the a band (r a = 120 ) and the /? 
band (T^ =60). Since the density of states (DOS) Npi 
at the Fermi level (_/V F 7 =0.53, iV Fa =0.10 and Np t 3 = 0 .37 
of the total DOS), the a band is neglected in this paper 
for simplicity. In the following we consider a two band 
model: the 7 band with _/V F 7 =0.53 and the /? band with 
7V F/ 3=0.47. 

In Sr 2 Ru 04 there are several outstanding unsolved is¬ 
sues: 

(1) Which bands is the major band for superconductivity, 
either the 7 band or the (3 band ? A^ > A 7 or A^ < A 7 ? 

(2) Can the PPE resolve the observed disparity be¬ 
tween the vortex lattice anisotropy^! 60 and the H c 2 
anisotropy!! 20 ? 


The first issue (1) has been extensively debated, in par¬ 
ticular in connection with the pairing mechanism of this 
material to stabilize the chiral p-wave stated—. We 
approach this issue from the different view point by an¬ 
alyzing the mixed state vortex states. The second issue 
( 2 ) has been investigated by u s 17 : 18 based on the single 
band picture. In this paper we revisit it based on a more 
realistic two band model, which is able to allow us to 
study the issue ( 1 ). 

We base our computations on the quasiclassical the¬ 
ory!! - —. The original single band theory is extended to 
various multiband cases, including MgB ; 22 : 23 , iron pnic- 
tides!! - — . The applicability of the quasiclassical the¬ 
ory!! - — is given in general by the condition fc F £ 1 with 
kp the Fermi wave number and £ the coherence length 
under the assumption that the normal state properties 
are described by a Fermi liquid theory. For Sr 2 Ru 04 
the three bands a , /3 and 7 are fc F =0.304, 0.622 and 
0.753 (A -1 ), respectively. The in-plane and c axis co¬ 
herence lengths £ are 660 and 33 (A), respectively!!. 
Thus for any combinations fc F £ 1 is well satisfied. 
We also notice that the lattice constants a=3.862A and 
c=12.722A are short enough compared with those coher¬ 
ence lengths, thus S^RuCU is a three dimensional nor¬ 
mal metal. Moreover, most physical quantities, including 
various transport coefficients and thermodynamic prop¬ 
erties can be consistently and coherently described by a 
Fermi liquid theory as explained in details by Mackenzie 
and Maeno!!. Therefore we can quite safely apply the 
quasiclassical theory to Sr 2 Ru 04 . 

By self-consistently solving microscopic quasiclassical 
Eilenberger equation with two bands, we calculate a vari¬ 
ety of the physical quantities relevant to available exper¬ 
iments, such as the form factors probed by SANS exper¬ 
iments*^!, magnetic torques!!, and the vortex lattice 
anisotropy Tvl which differs generally from the effective 
mass anisotropy mentioned above. 

The arrangement of this paper is as follows: After in¬ 
troducing Eilenberger theory extended to multibands!! 
with PPE in Sec. 2 which was done in our previous 
paper—, we construct a model system for considering 
Sr 2 Ru 04 by fixing the several model parameters in Sec. 
3. Here we compare the two scenarios, one is the /3 ma¬ 
jor and the other 7 in equal footing, finding that the 
former is better than the latter relative to the existing 
experiments. Then we come to the main theme of the 
present paper in Sec.4; computations and analyses of the 
form factors and vortex lattice anisotropies as a function 
of the angle D compared with the data of SANS exper¬ 
iments^*^. In Sec. 5 we examine the magnetic torques 
which are also measured recently by Kittaka et al!!. The 
final section 6 is devoted to discussions and conclusion. 
This paper is an extension of our previous work based on 
a single band mode l 18 : 29 : 30 and also closely related to our 
two band model calculation s 10 : 23 . 
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II. QUASICLASSICAL THEORY INCLUDING 
PPE FOR TWO BANDS 

We start with the free energy F in the quasiclassical 
thcoryi£~— extended to the two band case22r— , which is 
given by 

F = j ^{^^-^| J B(r)| 2 +E iJ A*(r)(U- 1 ) iJ A i (r) 

- nk B TN F o'Ei UJn \ <u]c '£j^p-(I(u} n ,k j ,r)) k } (1) 

A FO 

with Xn = 2^|7Vfo, N fo = T,jN Fj and 

I(u>n, kj,r ) = A(kj,r)f\uJn, kj,r)+A*(kj,r)f(cj„, kj,r) 

+ {9j - sgn(uj n )){^r(uj n +iiiBB + ^v F] ■ (V -i^-A))fj 
J j ^ @0 

+ ^ bB - + *X- A ))/|} (2) 

fJ 2 <Po 

The flux quantum <po = 2j£_. v F j is the Fermi velocity 
at kj of the band j. The Fermi surface average (• • • )k, is 
normalized within each band as (l)k = 1- Here we intro¬ 
duced the interaction matrix Vij with 2x2 for two bands 
where Vjj is the pairing interaction on the j band and 
Vij = Vji for i ^ j is the Cooper pair transfer between the 
i and j bands, gj = g(u n ,kj,r), fj = f(uj n ,kj,r), and 
/] = /'K , kj,r) are the quasiclassical Green’s func¬ 
tions for the j band. A j(kj, r) = Aj{r)(f>j(kj) is the pair 
potential and <f>j(kj ) describes the gap symmetry of the 
j band in reciprocal space, which allows us to choose the 
gap form depending upon each band as will be done in 
the following. The vector potential A(r) and the internal 
field B(r) are related to B(r) = V x A = B + b(r) with 
B uniform field. 

By following the same procedure by Eilenbergei— , the 
functional derivatives with respect to fj = f(co n ,kj,r ), 
and f) = / f (w n ,kj,r) yield the so-called Eilenberger 
equation extended to the two band case: 

{uj n + ig,B{r) + vj ■ [V + iA(r)]} fj = A j (k j ,r)g j , 

{u n + igB(r ) - vj • [V - U4(r)]} f) = A* j (k j ,r)g j . 

( 3 ) 

This form is understandable because the fourth term 
in the free energy Eq. ©, which includes the Green’s 
functions, is separable in the band index, thus the re¬ 
sultant equation of the functional derivative should be 
separable for each band. The stationary conditions of 
Eq. © with respect to the functionals A *(r) and the 
vector potential A(r) give rise to a complete set of the 
self-consistent equations extended to the two band case, 


which are given below; Eqs. © and ©. This complete 
set of the self-consistent equations coincides and is con¬ 
sistent with those obtained previousl y 22 ! 23 . 

The electronic state is calculated by solving the Eilen¬ 
berger equation Eq. © in the clean limit^I, including 
the Pauli paramagnetic effect (PPE) due to the Zeeman 
term /rH(r) 32 , where g, = /iBBo/nkETc is a renormalized 
Bohr magneton related to the so-called Maki parame¬ 
ter a m = 1.76/r. The quasiclassical Green’s functions 
g{kj, r, uj n + ifj,B), f(kj,r,u) n + ijiB), and / t (fc J ,r,w n -1- 
igB) with the band index j depend on the direction of 
the Fermi momentum kj for each band, the center-of- 
mass coordinate r for the Cooper pair, and Matsubara 
frequency uj n = (2n+l)7rfcBT with n€Z. They are calcu¬ 
lated in a unit cell of the triangle vortex lattice. 

The unit of Fermi velocity ufo is defined by -/Vfo^fo = 
Afi^pi T A^f 2 ^f 2 ) where the density of states (DOS) in 
the normal state at each Fermi surface is defined by 
TVfo = Api + A r p 2 - Throughout this paper, tempera¬ 
tures, energies, lengths, and magnetic fields are, respec¬ 
tively, measured in units of the transition temperature 
T c , nkBT c , £o = /w F o/27rfc B 7c, and B Q = 4> 0 /2it£ ! l. We cal¬ 
culate the spatial structure of g in a fully self-consistent 
way. 

The pairing potential Aj(r) is calculated by the gap 
equation 

Aj( r ) = t J2 E {Ui + f}*)Mh)) k , 

0<UJn<UJ c £= 1,2 

( 4 ) 

which is coupled via the interaction matrix V. We use 
the energy cutoff w c = 20fcBT c . The vector potential is 
also self-consistently determined by 

V x V x A = V x Aip ara — —j N F j(vjlm[gj]) k ^, 

l^nl<“ci=l> 2 

( 5 ) 

which includes the contribution of the paramagnetic mo¬ 
ment iV^Zp ara = (0, 0, -A^/para) with 

^7 para = M 0 | —g —jy E E ^ ) ■ 

\ ^ |Wn|<WcJ=l,2 / 

( 6 ) 

Here B is the averaged flux density mentioned above, the 
normal state paramagnetic moment Mq = (g,/ k) 2 B, and 
k = Bo/TrkBT c y/8TrN F o. The Ginzburg-Landau (GL) pa¬ 
rameter kql is the ratio of the penetration depth to co¬ 
herence length for B || c. Using Doria-Gubernatis-Rainer 
scaling^, we obtain the relation^ of B and the external 
field H. The total magnetization M tota \ = B — H includ¬ 
ing both the diamagnetic and the paramagnetic contri¬ 
butions is derived. 

We solve Eq. © with ito n —>• E + irj for the electronic 
state. The local density of states (LDOS) is given by 
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Nj(r, E) = Nj^(r, E)+Nj > i(r, E ) where 
Ay c t( r, E) = Af’j (Re ^g{kj , v , -t -/<t / f/i)| iuj n —t&Hr]]') k ^ ■ 

V) 

with er = 1 (—1) for up (down) spin component. We 
typically use the smearing factor 77 = 0.01. The DOS is 
obtained by the spatial average of the LDOS as N(E) = 
Yhj Nj{E) = J2j(Nj^( r , E)+Nj^(r, E)) r . 

We consider a simplified model of a two-band sys¬ 
tem with a larger superconducting gap band (band- 
1) and a smaller gap band (band- 2 ). As a model 
of the Fermi surfaces, we use two quasi-two dimen¬ 
sional Fermi surfaces with rippled cylinder shapes. The 
Fermi velocity is assumed to be Vj = {vj t a,Vj t b,Vj tC ) oc 
(coscf), sin(j>, Vj }Z smkj tC ) at kj = {kj t a:kj t b,kj tC ) oc 
(kj cos <p. kj sin <p. Vj tC ) on the Fermi surfaces^. We con¬ 
sider a case Vj, z = 1 /Tj, to produce large anisotropy ratio 
of the coherence lengths, 



FIG. 2: (Color online) The specific heat data (open sym¬ 
bols)^ analysis by the /? major scenario (bold line) with 
A^(0)/A 7 (0) = 2.5 where the line node (full) gap is on the 
7 (p) band and by the 7 major scenario (dotted line) with 
A 7 (0)/A^(0) = 1.7 where the both bands contain the line 
nodes. For comparison, the standard BCS case with the full 
gap is also shown (thin line). 


r l = &>/&«. ^ <4c>kf /<4b>kf 


( 8 ) 


with j = 1,2 where (• • ■ } kj indicates an average over 
the Fermi surface on each band. The magnetic field ori¬ 
entation is tilted by 9 = 90° — 11 from the c axis to¬ 
wards the ab plane. Since we set z axis to the vortex 
line direction, the coordinate r = (x, y , z) for the vortex 
structure is related to the crystal coordinate (a, b , c) as 
(x, y, z) = (a, b cos 9 + csm9,ccos9 — bsin9). 

We set unit vectors of the vortex lattice as 


Ul = c(a/2, —\/3/2), u 2 = c(a/2, y/3/2) (9) 

with c 2 = 2</> 0 / (V3aB) and the vortex lattice anisotropy 
is defined by Tvl = a/y/P. The anisotropic ratio 

r #) = tj,v/ij,x ~ ( v ly)H 3 2 /(^,x>kf 1 which comes from 

the Fermi velocity anisotropy, 

r#) = , 1 „ = (10) 

J cos 2 9 + 2 sin 2 9 


with j = 1 and 2. 

To discuss R-dependence of the internal field distribu¬ 
tion B(r) = V x A, we consider flux line lattice (FLL) 
form factor F(q h>k ) = (F X (M)> which is 

obtained by Fourier transformation of the internal field 
distribution as B(r) = J2 h ,k exp^q^*, • r) with 

the wave vector = ftqi + fcq 2 . h and k are inte¬ 
gers. The unit vectors in reciprocal space are given by 
qi = (27r/c)(l/a,-l/v / 3) and q 2 = (2iv/c)(l/a, l/y/S). 
The ^-component | F z ( k ,k) | 2 from B z (v) gives the intensity 
of conventional non-spinflip SANS. The transverse com¬ 
ponent, |F tr(h;k) | 2 = |F x(M) | 2 + \F y(htk) \ 2 , is accessible 
by spin-flip SANS experiments ! 12 ' 35 


III. MODEL CONSTRUCTION AND PHYSICS 

FOR H || c 


A. Specific heat and H c 2 anisotropy ratio 


In order to determine the model parameters appropri¬ 
ate for Sr 2 RuC> 4 , we start out to fix the gap magnitudes 
A 7 (T) and Ap(T) and their nodal structures. We first 
analyze the electronic specific heat data C{T)/T at zero 
fields by solving self-consistently the Eilenberger equa¬ 
tion Eq. m for the uniform system without using the 
phenomenological so-called a model^. Here we have as¬ 
sumed that the minor component induced only by the 
Cooper pair transfer coupling Vi 2 . The direct attrac¬ 
tive coupling among the minor component is vanishing; 

F 22 = 0. 

As seen from Fig [2] the C(T)/T data is equally well 
explained either by 

(A) the 7 scenario with A 7 (0)/A l g(0) = 1.7 at T = 0 
where both bands have line nodes as coincided with other 
authors^, 

(B) the /3 scenario with A l g(0)/A 7 (0) = 2.5 where the j3 
band has a full gap and the 7 band line nodes. 

It is apparent that the linear T behavior of C(T)/T at 
lower T indicates that the nodal gap is necessary some¬ 
where; for the 7 scenario in both bands and for the the 
P scenario only in the 7 band. This nodal structure dif¬ 
ference between them is decisive for the two scenarios as 
will be seen shortly. 

According to Kittaka et al2£, the upper critical field 
ratio H C 2 ,ab(T)/H c2t c(T) for the two directions (H || ab 
and H || c) is T-dependent, implying that PPE becomes 
stronger as the field applied to the ab plane increases. 
Here we calculate Lf c 2 )a b(T)/7L c 2 iC (T) for both scenarios 
and depict the results in Fig. [3] Since near T c the in¬ 
trinsic effective mass anisotropy governs its tending limit 
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( a ) F z (io/^c2, c (b) F z (ii)/W C 2 j c 




FIG. 3: (Color online) The upper critical field ratio 
H C 2 ,ab(T)/ H C 2 ,c(T) for the two directions (H || ab and H || c). 
The two cases for the /3 (filled triangles) and 7 scenarios (filled 
squares) are compared with the experimental data (open sym¬ 
bols)^. This shows that the P scenario is superior to the 7 
scenario. 



FIG. 4: (Color online) The comparison of the experimental 
data (open symbols)^ of the specific heat at T=60mK as a 
function of J7(|| c) with the theoretical results (filled symbols) 
for the P and 7 scenarios (T = 0.1T C ). This shows too much 
low energy excitations released at lower fields in the 7 scenario 
because the nodal gap exists in the minor band. 

(T ->• T c ), namely, H c2%ab /H c%c -)■ 180 (60) for the 7 (/?) 
scenario. The experimental data££ shown support the /3 
scenario within the experimental accuracy where there 
is no indication of the ratio with tending to 180. This 
is one of the most clear signatures for the /3 major sce¬ 
nario. By adjusting the /1 parameters for both scenarios 
the best fittings are accomplished in /r=0.04 (0.02) for 
the /? (7) scenario. From now on we use those values in 
the following calculations. 

In order to further distinguish between the two sce¬ 
narios, we take up the experimental data£& of the field 
dependence of C(H)/T for H || c at the available low¬ 
est temperature T=60 mK, which best mimics the zero- 
energy density of states N(H) at T = 0 in the theoreti¬ 
cal calculation. In Fig. [4] we compare the experimental 
data with the theoretical values as a function of the field 
applied parallel to the c axis. As seen from it, the 7 
scenario overestimates the data while the /3 scenario un¬ 
derestimates it. If taking into the finite T effect in the 


FIG. 5: (Color online) The longitudinal form factors F z (io) 
and F z (n) as a function of J7(|| c). The experimental data 
(open symbols are compared with the two scenarios. The 
theoretical results (filled symbols) are multiplied by a factor 
1.7 (3.5) for the p ( 7 ) case (T = 0.1T C ). This shows that the 
P scenario is superior to the 7 scenario. 


experimental data, the theoretical curves should move up 
when considering thermal excitations, resulting in fur¬ 
ther departure of the 7-major curve while the /3-major 
curve comes closer to the data. We also note from Fig. 
U that too much low energy excitations are released at 
lower fields in the 7 scenario because the nodal gap ex¬ 
ists in the minor band, causing the overestimate. This 
fact, which has been unnoticed so far is quite fatal for 
the 7 scenario. 

Thus it is clear from those two criteria that the /3 sce¬ 
nario is far better than the 7 scenario. The best fitting 
also yields the kql value for H || c, that is, kgl=2.7 
whose value is used in the following computations. But 
we should keep in mind that this low kgl value delicately 
depends on the particular sample used. This quantity is 
known to be sample dependent. 


B. Form factors for H || c 

The form factors of the longitudinal components -F'z(io) 
and iTj/n) which are measured by SANS experiments^ 
for H || c are compared with the two scenarios whose 
magnitudes are multiplied by a factor of 1.7 for the /3 sce¬ 
nario and 3.5 for the 7 scenario. This is partly because 
the actual kgl= 2.7 determined by C(H)/T fitting previ¬ 
ously might be different from the SANS experiment (In 
fact 37 C 2,c = 58mT and kgl=2.0 differ from H c 2 iC =75mT 
and kgl=2.3 in the best samplesii). Since the form fac¬ 
tor magnitudes are sensitive to the kgl value (oc Kq^), 
it is permissible to adjust it to fit with data. As seen 
from Fig. [5] the field dependences of Fb(io) and ^b(n) 
are far better explained by the /3 scenario than by the 7 
scenario. 

In Fig. [6] we decompose the longitudinal form factors 
F z ( 10) and -Fqu) into the two contributions of the /3 band 
and 7 band in the case of the /3 scenario. It is seen that 
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Fig. 1(c) at 12 = 0° respectively. Note that the crossing 
point A of the two orbital limited H c 2’s in Fig. 1(a) is 
located at around Ha ~ 1°. 


IV. FORM FACTORS AND VORTEX LATTICE 
ANISOTROPY 

A. Form factors 


FIG. 6: (Color online) The longitudinal form factors F 2 (io) 
and F z (n) are decomposed into the the /3 band contribution 
and 7 band contribution in the case of the /3 scenario (T = 
O.lTc). 


the minor 7 contribution amounts to ~ 10% of the total 
at this temperature T = 0.1T C for both F z ( 10) and FGn). 
Thus in order to understand the field dependence of the 
form factors, the multiband effect is essential, which fur¬ 
ther becomes clear later. Note in passing as mentioned 
in Introduction on MgB2 Cubitt et a£ discover the ad¬ 
ditional contribution of the minor 7r band to the main a 
band contribution at lower H. This general trend here 
supports their discovery (see Fig. 2 in Ref. [§!). Note 
that the relative weight of the main /? contribution and 
minor 7 contribution depends on temperature, field and 
field orientation 12. Generally as T and H decrease, the 
minor 7 contribution increases because the two order pa¬ 
rameters are more competitive there. 


C. Phase diagrams in H vs 12 

Since we have determined the gap ratios for each sce¬ 
nario, it is possible to establish the phase diagrams on H 
vs 12 plane. For H \\ c the two orbital limited H l c2 c ratio 
(* = /?, 7) with the obvious notations is written as 




H 1 

n c2 ,c 


£2 

^7,c 

— 

“SySjC 



/ 

W,C 


2 


( 11 ) 


It is knownlA that the Fermi velocity ratio v~ tC /vp tC = 
0.5. Thus H^ JHl 2 iC = 2.5 2 x 0.5 2 ~ 1.6 for the /3 
scenario while H^ 2c /H^! 2c = 1.7 -2 x 0.5 2 ~ 0.1 for the 7 
scenario, those corresponding to Fig. 1(a) and Fig. 1(c) 
at 12 = 90° respectively. On the other hand, at 12 = 0° 


ff'/ 

n c2 ,ab 

H 13 

n c2 ,ab 


T 1 H c2,c 

r P H c2, C 


= 3 X 


H 1 

I1 c2,C 


H. 


c2,c 


( 12 ) 


This yields H 3 2 ab /H^ 2 ab = 3/1.6 ~ 1.9 for the /3 scenario 
while H^ 2 ab /H^ 2 ab = 3/0.1 ~ 30 for the 7 scenario, thus 
giving rise to the situations depicted in Fig. 1(a) and 


The form factors FF are a sensitive and useful probe 
measured by SANS in order to detect the Held distribu¬ 
tion in the vortex state. The SANS experiment s 12 ’ 27 on 
Sr2Ru04 are performed. They find the transverse com¬ 
ponents of FF as a function of 11 near the ab plane. Here 
we obtain FF by evaluating the field distribution via a 
self-consistent solution of Eq. ©. 

It is seen from Figs. [7] (a) and (b) that the angle 12 
dependences of the transverse component of FF for two 
fields B = 2.0 (high field) and B = 0.5 (low field) ex¬ 
hibit a different characteristic; The high field data show 
a simple monotonic decrease after taking a maximum 
towards 12 = 12 c at which the superconducting-normal 
state transition takes place. The monotonic transverse 
FF curve is similar to that of the single band case shown 
in previous paper s 18 ’ 30 . This high field scan (B = 2.0) 
corresponds to the horizontal scanning path in Fig©b) 
where only the Tp region is sensed. The maximum posi¬ 
tion 12^ x = 0.9° coincides with those of the single band 
result^ with T = 60. The overall features of the experi¬ 
mental data for H = 0.5T and 0.7T are well reproduced 
as seen from Fig. [7] (a) except for a few data points at 
higher angles. 

On the other hand, for the low field result (B = 0.5) 
shown in FigJT] (b), it is seen that the theoretical FF 
curve coincides with the maximum angle of the experi¬ 
mental data at H = 0.15T and 0.25T. However, it devi¬ 
ates from those data after that. The experimental data 
tend to vanish at round 12 c ~ 6°, which is by far from the 
known 12 c ~ 30°(18°) at H = 0.15 (0.25)T, which should 
be ultimate vanishing angle for FF amplitude. We will 
discuss this discrepancy shortly. Here we just point out 
that those low field data correspond to the low field scan¬ 
ning paths in FigOJb) where the crossover of the regions 
—> T 7 —> r,3 is sensing with increasing 12. 

We decompose the FF contributions from the major 
13 band and the minor 7 band for the transverse compo¬ 
nents Ft r (io) ( a ) an d -Ftr(H) (b) and the longitudinal com¬ 
ponents F z ( 10 ) (c) and F z ( 10 ) (d) as shown in Fig.[8]for the 
high field and Fig.[!j]for the low field. The relative weight 
of the minor component in this transverse FF is around 
10% of the total FF and their peak contribution coincides 
with the maximum position f2 max for both B = 2.0 (high 
field) and B = 0.5 (low field) data. As for the longitu¬ 
dinal components of F z ( 10) and F z ( 10 ), it is reasonable to 
see that the maximum angle 12 max is shifted to a higher 
angle than F’ tr ( 1 i). This is because when i7 C 2(S2 c ) is ap¬ 
proached in which the “effective” field virtually increases 
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FIG. 7: (Color online) The transverse form factors Ftr(ii) of 
the SANS data (open symbols ) 12 ' 27 at 22=0.7T and 27=0.5T 
are compared with the theoretical results (filled symbols) for 
B = 2.0 (a) and the data at H=0. 25T and 0.15T with B = 0.5 
(T = 0.527) (b). 
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FIG. 8: (Color online) The decomposition of the transverse 
components of the form factors Ftr(io) ( a ) an d Ftr(ii) (b) at 
high field B = 2.0 (T = 0.527). The decomposition of the 
longitudinal components of the form factors 27,(10) ( c ) an d 
2^(11) (d) at high field B = 2.0 (T = 0.5T C ). 


because the decreasing A (ST) means decreasing H p and 
hence enhances PPE, thus pushing up the longitudinal 
FF \F Z \. 

In order to clearly see the different angle dependences 
Ftr(n) (12) for the two cases and to understand the dis¬ 
crepancy of the low held data mentioned above shown in 
Fig. [3(b), we replot those theoretical results normalized 
by the own f2 C in Fig. llOl and compare those with the cor¬ 



FIG. 9: (Color online) The decomposition of the transverse 
components of the form factors 27 r (io) ( a ) an d 2 ? tr (n) (b) at 
low field B = 0.5 (T = 0.527). The decomposition of the 
longitudinal components of the form factors F x iiq\ (c) and 
27,(11) (d) at low field B = 0.5 (T = 0.527). 


responding single band results^ where we have adjusted 
the vertical scale so that the slopes of these curves near 
f l c coincide with each other. It is now seen clearly that 

(1) The high field result B = 2.0 belongs to the single 
band universality curve. This is because the high field 
scanning path is sensing only the T p region in Fig. Q] (b). 
This is virtually same as in the single band case. Namely 
after taking the maximum the FF curve simply goes to 
vanish at f2 C . 

(2) The low field result B = 0.5 behaves differently from 
those single band universality curves and exhibits a “bi- 
mo dal” 12 dependence where the FF peak is additionally 
enhanced. Thus the curve just after taking the maximum 
tends to vanish earlier than at Q c which is the ultimate 
vanishing angle. 

Let us come back to understand the low field FF data 
shown in Fig. [3(b). The bimodal F tr (n)(Q) structure at 
the low field shown in Fig. [3 (b) where the slope just af¬ 
ter the maximum differs from the slope at higher angles 
as increases. This crossover angle I2 cross ~ 4° — 5°. 
This low field scan corresponds to the low field scan¬ 
ning path in Fig. 1(b) where the crossover of the regions 
—> T 7 —> is sensing with increasing f l. In partic¬ 

ular, the central peak II region (1.0° < fi < 4.0°) in the 
bimodal structure corresponds to the middle T 7 region in 
this crossover. This theoretical curve is compared with 
the SANS data^l at H = 0.25T in Fig[3 (b). Although 
































FIG. 10: (Color online) The comparison of the transverse 
form factors F tl .m) with the single band results (open sym¬ 
bols)^ for S=1.5, 3.0, and 4.5 and the present multiband 
results (filled symbols) depicted in Fig. [7] for 5=2.0 and 0.5. 
The vertical scale is adjusted so that the slopes of these curves 
near S2 C coincide with each other. It is clear that the low field 
theoretical result (5 = 0.5) behaves differently, exhibiting a 
bimodal structure. 


the present experimental data do not exhibit this bimodal 
FF structure, the existing data are understood as coming 
only from the middle T 7 region among the bimodal struc¬ 
ture because the extrapolated vanishing angle fl c ~ 7°. 
This critical angle is too small compared with f l c = 18° 
which must be the ultimate vanishing angle of the FF 
intensity. Thus we interpret this so that the existing FF 
data are coming for the peak T 7 region of the bimodal 
FF distribution. This interpretation is strengthened later 
when analyzing the vortex lattice anisotropy Tvl- 


B. Vortex lattice anisotropy Tvl 

Figure [TT] shows the results of the vortex lattice 
anisotropy Tvl compared with the SANS experiments^ 
where for a given B and ft the self-consistent solutions 
of Eq. Q are optimized for varying Tvl as shown in 
Fig. [TT] to seek the free energy minimum where the free 
energy form is given by Eq. ©■ In general the opti¬ 
mized Tvl is that given neither by the effective mass 
anisotropies F^fi) with r / 3=60 nor F 7 (fi) with F 7 =180 
alone given by Eq. m- This is obvious because those 
two anisotropies are coupled and competed by the multi¬ 
band effect. 

The higher field results in Fig. [11] (a) follow rather 
well those given by the above formula of Eq. (fTOl) for 
Tp(n) for the higher angles > 3°. We note that at 
ft = 0 ° where Fvl( 0°)=70 corresponds precisely to the 
anisotropy H^ 2 (0°) / H^ 2 (90°) for T = 0.5 T c . Between 
the angles 1° < ft < 3° the theoretical results devi¬ 
ate upwards, that is, Fvl(^) > r^fi). The high field 



FIG. 11: (Color online) The comparison of the experimen¬ 
tal data (open symbols) 12 ' 27 for Fvl(B) for H = 0.5T with 
the high field results of 5 = 2.0 (filled symbols) (a) and 
H = 0.25T with the low field results of 5 = 0.5 (filled sym¬ 
bols) (b) (T = 0.5T C ). The curves are drawn by the effective 
mass formula Eq. (1101) with T = 60 (dotted line) and F = oo 
(dashed line). 


scan in Fig. 1(b) is barely touching the T 7 region whose 
anisotropy can be certainly larger than r j a=60 because 

r 7 =180. 

The low field theoretical results shown in Fig. [TT] (b) 
follow the F j a(n) curve for 0° < ft < 1° and remark¬ 
ably exceed the T —> oo curve in 1° < ft < 5°. This 
window in ft corresponding to the T 7 region appears be¬ 
cause there the F 7 =180 anisotropy further modifies the 
Tvl value upwards, simultaneously enhancing the trans¬ 
verse FF. We point out a theoretical fact that as ryL(n) 
increases, F tr (Q) becomes larger. We anticipate that the 
results at lower temperatures than the present one at 
T = 0.5 T c would improve the quantitative fittings of TyL 
and simultaneously F tI (fl). Note that ryL(fI = 0°)=52 
corresponds to the single band anisotropy for the /3 band 
at lower fields^. The H = 0.25T data nicely fit our 
theoretical result. The important point here is that not 
only the experimental data exceed the T —> oo line, but 
also there exists a wider window 1° < ft < 5° where the 
experimental data deviate from the single band F^(fl) 
curve. 

We demonstrate in Fig. [TT] that the optimal Fvl is 
determined careful enough, which can be larger than the 
corresponding T —> oo case in Eq. (HDD- For example, for 
our 5=0.5 case at ft = 2° our Tvl= 32 while Tvl= 29 for 
r = oo and FyL=28 for T^ = 60. Thus it can be said 
that the multiband effect helps enhancing TyL beyond 
r = oo. 


C. Order parameters, free energy and 
magnetization 

We illustrate the ft dependences of several physical 
quantities of interest which are the basis of the form fac¬ 
tors and magnetic torque calculations as shown shortly. 
The free energy F(ft) is shown in Fig. [T3] from which we 
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FIG. 12: (Color online) The free energy curves for various 
angles Cl as a function of Fvl for B= 2.0 (filled symbols) and 
13=0.5 (open symbols) in order to emphasize that our compu¬ 
tations are performed accurately enough to resolve the subtle 
differences of the Fvl values. 

will evaluate the magnetic torques. Since the transition 
at Cl c is of second order at those fields and temperature 
(T=0.5T C ), F(Cl) becomes zero smoothly at the transi¬ 
tion point Cl c . 

As seen from Fig. Ildl the two order parameters; the 
major A l g(f2) and minor A 7 (f2) start decreasing from 
Cl = 0° towards Cl c as a function of Cl. The two 
curves change smoothly in parallel because the minor 
component A 7 (f2) is induced by the major one A l a(f2) 
through the Cooper pair transfer V 12 in the absence of 
V 2 2 - In this case we expect no hidden first order transi¬ 
tion phenomena^ where there is an abrupt change of the 
two order parameters inside the superconducting state. 

In Fig. [15] the total paramagnetic susceptibility xt(^) 
and decomposed Xfi(Cl) and X 7 (C) are displayed as a 
function of Cl. Since we assume the density of states 
N-pry = 43% and Npp = 57% of the total DOS, the cor¬ 
responding paramagnetic values X 7 (D) > Xp(Cl) as ex ~ 
pected. We also point out that the base paramagnetic 
moment values at Cl = 0° is large because the calcula¬ 
tions are done at rather high temperature T = 0.5T C . 
Needless to say, at T—► 0 xt should vanish at Cl = 0°. As 
Cl increases y t (0) becomes larger because the system is 
approaching the transition point at Cl c where the normal 
value xn must be recovered, that is, Xt(flc) = Xn- 

V. MAGNETIC TORQUES 



FIG. 13: (Color online) The angle dependences of the free 
energy for 13=2.0 (filled symbols) and 13=0.5 (open symbols), 
showing a smooth change of the second order transition at 
where Cl c = 7.78° and Cl c = 33.96° respectively. 


A/u%r c 



FIG. 14: (Color online) The angle dependences of the order 
parameter amplitudes Ap(Cl) and A 7 (fi) at the center of the 
vortex unit cell for 13=2.0 (filled symbols) and 13=0.5 (open 
symbols). 



Magnetic torque t(CI) defined by t(CI) = dF(Cl)/dCl 

provides several important information on a uniaxial FIG- 15: (Color online) The angle dependences of the para- 
anisotropic superconductor. We can know the intrinsic magnetic susceptibility Xt(D) for 13=2.0 (filled symbols) and 
anisotropy of a system by the peak position of the torque S=0.5(open symbols), which are decomposed into Xfl (Cl) and 
curve t(CI). This can be easily performed by using a phe- '' 
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(a) x(n) [a.u.] (b) x(n) [a.u.] 



FIG. 16: (Color online) The angle dependences of the torques 
(filled symbols) for 13=2.0 together with the experimental 
data (open symbols)^ of H=0. 5T (a) and 13=0.5 with the 
data of H=0.2T (b) 


nomenological theory based on London theory^ for sin¬ 
gle band superconductors. In our previous paper s 18 ’ 30 we 
examine the applicability of this approach and propose 
a modification to this. Extending this single band Eilen- 
berger calculations, here we show the results of the torque 
curves r(ff) for the present two band model and use those 
to analyze the data on Si^RuO-r where the torque curves 
are measured recently^. 

We first display our results for the torque curves in 
Fig. [TCI where our results for B = 2.0 and B = 0.5 are 
compared with the experimental data for H= 0.5T (for¬ 
mer) and H= 0.2T (latter). It is seen from Fig. [16] (a) at 
the high field B = 2.0 that the fitting is done well, such as 
the peak position and the Yl c value. Since the high field 
result is sensing only on the Yp region, it is reasonable 
that the theoretical curve nicely explains the experimen¬ 
tal data at F=0.5T and the peak position = 1-5° 

coincides with the single band result^. This peak posi¬ 
tion Hp° r X e also coincides with HpJ ak of the form factor 

(B = 2.0). Thus the four values Hp° r X e and ^pf a k both 
for theory and experiment are coinciding with each other, 
leading us to firmly conclude that in high fields the sys¬ 
tem is virtually in the single band-like T p region. 

On the other hand, the low field result (B = 0.5) shown 
in Fig. [TC] (b) explains the experimental torque curves 
for L/=0.2T, but the peak position differs slightly from 
each others. We also notice here that the theoretical 
0 U V^ e F ak and ^peak (~ 4°) agrees with the experi¬ 
mental data as shown before. 

It is interesting to compare our torque curves with 
those for the single band case which are displayed in 
Fig-HZl It is seen that the two torque curves nicely corre¬ 
spond to the single band curves, which is in contrast with 
the FF case shown in Fig. [TO] where the low field curve 
markedly deviates from the single band case. This means 
that the magnetic torque is a rather insensitive probe to 
see the subtle, but important multiband effect. In other 
words, the form factor measurement is sensitive enough 
to distinguish the multiband effect from the single band 
effect. This is because the torques r(fl) = dF(Cl)/dCl 



FIG. 17: (Color online) The torque curves as a function of 
Cl for 13=2.0 and i3=0.5 at T = 0.5T C . Those are compared 
with the single band results^ for i3=1.5, 3.0, and 4.5. The 
vertical scale is adjusted so that the slopes of these curves 
near Cl c coincide with each other. 
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FIG. 18: (Color online) The comparison of the Fvl with the 
single band cases Tvl(60) and Fvl(180) and also with the 
H c 2 anisotropy H C 2 (Cl) / H C 2 , C - (a) High field case 2.0 where 
those three data Tvl, Fvl(60) and H C 2 /H C 2 , C coincide with 
each other, and (b) low field case 0.5 where the lower (high) 
angle region corresponds to the /3 ( 7 ) band anisotropy. The 
vertical dashed line at around Cl = 1.5° denotes its boundary. 
The dotted (dashed) curves indicates the angle dependence of 
the effective mass formula Eq. (1101 for T = 60 (r = 00 ). 


comes from the total free energy F(Cl) while the FF is 
probing the particular Fourier component of the spatially 
modulated magnetic field in the mixed state selectively. 
Thus it is natural to expect that the FF measurement is 
more sensitive than the torque measurement in picking 
up the multiband effect. 


VI. DISCUSSIONS 

A. single band vs multiband 

In order to describe the transverse components of FF 
and Tvl data32 near the ab plane which are the main 
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themes of the paper, we have done the calculations based 
on both scenarios, first focusing on the H || c physics 
to fix the model parameters. Here we critically exam¬ 
ine those scenarios comparatively. It is obvious that the 
multiband scenario superior than the single band sce¬ 
nario^ because the former includes the latter as a lim¬ 
iting case. A question is how effectively the single band 
scenario can describe those data or how the multiband 
description is inevitable for the Sr2Ru04 physics. As al¬ 
ready demonstrated and also shown in Fig. [T8l the over¬ 
all features of those data can be reproduced by the sin¬ 
gle band model with T j a=60 in the higher fields above 
H > 0.5T in the H vs Q plane where there is a lit¬ 
tle trace for the existence of the 7 band with T 7 =180 
(see Fig. UsT a)). In this sense the single band picture is 
enough for this region. In contrast, however, it is neces¬ 
sary to retain the 7 band contribution in addition to the 
major /3 band contribution in the low field region below 
0.5T (see Fig. ITST bY). The former contribution is hid¬ 
den and not explicit where the crossover occurs around 
f l = 1.5° as indicated in Fig. ITBI bb It is necessary and 
indispensable to take it into account for explaining the 
FF and Tvl data—. This is a part of the reasons why 
our single band theory^ is successful to understand the 
mixed state properties of the present strongly uniaxial 
anisotropic superconductor. Moreover, depending upon 
the physical quantities of interest, the multiband effect is 
not so apparent even in the low field region. Namely, the 
torque curves are quite insensitive to the presence of the 
minor band and all the theoretical torque curves collapse 
into a universal curve by appropriate scaling as discussed 
before (see Fig. ED. 


B. ft main scenario vs 7 main scenario 


The question on which band the gap is large is much 
discussed in connection with the pairing mechanism to 
stabilize the chiral p-wave state by many authorsHvIli. 
Here we take a different view to consider this question 
through the analysis of the FF, Tvl and torques. If the 7 
band is major, this survives in the high field region of the 
H vs Q plane over the minor /3 band. Then the SANS ex¬ 
periment should detect it, manifesting itself in Fvl, that 
is, Tyl —> 180 as Yl —> 0. However, both high fields and 
low filed data shown in Fig. fill indicate Tvl ~ 60 or so, 
which is direct evidence that the /3 band with Yp = 60 
is major. This conclusion is also supported by the H c 2 
ratio H c2 ,ab(T)/H C 2 t c{T) 60 when T —> T , ^ because 
near T c this ratio directly reflects the intrinsic Fermi ve¬ 
locity anisotropy^, or the coherent length anisotropy. 
Thus those experimental data obviously reveal that the 
/3 band is major with certainty. 


C. Cooper pair tunneling V 12 vs direct attraction 
V 22 : Hidden first order 

In general the two band contains the three pairing pa¬ 
rameters Hu (V22) is the attractive interaction for the 
major (minor) band, and V12 = V21 is the Cooper pair 
tunneling term or proximity induced term. Here we set 
1/22=0 in this paper because there is no or little indication 
for V22 7^ 0, which causes the so-called hidden first order 
phenomena—. Namely at H* inside H c 2 certain physi¬ 
cal quantities exhibit a sudden change as a function of 
H , such as the Sonnnerfeld coefficient 7 (H) or the mag¬ 
netization curve as observed in CeCr^Si^, UBei^ and 
KFe2As2^- This means that V22 ~ 0 in S^RuCH and the 
minor 7 band pairing is exclusively induced by the major 
7 band through V12. 

D. Expected low T behaviors 

Because of several technical reasons our main com¬ 
putations have been done at relatively high T, namely 
T = 0.5T c , which are nevertheless successful in capturing 
the characteristic features in the FF, Tvl and torques in 
a qualitative level. To understand those data in a quan¬ 
titative level it is necessary to go into lower T which is 
demanding computationally. Here we anticipate possible 
outcomes if performing it. As shown in Fig. [7] we have 
mentioned that the low filed FF (13=0.5) differs from that 
in 13=2.0, which belongs to the single band universality 
class. The FF in low angle Q is enormously enhanced 
by the assistance from the minor band which amounts 
almost 10% even at T = 0.5T C . In low T calculation this 
enhancement of the FF should increase and the FF an¬ 
gle dependence becomes more similar to the data shown 
in Fig. Q] (b). This expectation is reasonable because at 
lower T the induced paramagnetic moments are confined 
in the vortex core, making more contrast in the spatial 
field distribution and thus enhancing the FF amplitude. 
See also Fig. fl5l where the paramagnetic susceptibility or 
paramagnetic moments of the minor 7 band are almost 
exhausted at H=0 to its normal value, meaning that is 
not confined in the core. Thus, the paramagnetic mo¬ 
ments are spreading out the whole space uniformly. 

E. Other multiband superconductors with PPE 

We can deepen our understandings of the present ma¬ 
terial by comparing similar multiband superconductors, 
CeCh^Si^, UBeigi and KFe2As2^- The first two are 
heavy fermion materials known to have multiband with 
all full gaps belonging to the spin singlet category where 
H c 2 is strongly suppressed and the hidden first order 
like anomalies at H* exist. Thus as mentioned V22 is 
indispensable for those systems. According to the re¬ 
cent SANS experiment^ the typical multiband Fe pnic- 
tide KFe2As2 is similar to Sr2Ru04 in the point that 
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the vortex anisotropy Tvl observed differs from the H c2 
anisotropy. This can be also explained by PPE. Thus 
the present theoretical framework, which is quite general 
might be able to explain those systems. 


F. Gap structure: vertical line node vs horizontal 
line node 

We have assigned the nodal structure that the domi¬ 
nant /? band is a full gap while the minor 7 is vertical 
line nodes. In order to reproduce the square vortex lat¬ 
tice oriented to the (110) direction observed for H || c 
which dominates the whole space measured in the H vs 
T phased. The gap (or near) nodes should be in this di¬ 
rection in reciprocal spaced, implying a d x 2_ y 2 like nodal 
structure, contrary to the claim by Deguchi et al2£ who 
measure the specific heat by rotating the field direction 
and see the four-fold oscillation patterns whose minima 
are located for the (100) direction. They conclude that 
the 7 band has a d xy like nodal structure. However, this 
assignment is difficult to explain the square vortex lattice 
orientated along the (110) direction for H || c. Even if we 
take into account the in-plane Fermi velocity anisotropy 
in the /? and a bands whose minima are oriented along the 
( 100 ) direction^, thus preferring the square lattice ori¬ 
ented along the ( 100 ) direction as discussed previously— 
Thus the best way to avoid this difficulty is to simply con¬ 
sider that the 7 band nodal structure is d x 2_ y 2 like. The 
specific heat oscillation experiment done above 120mK 
is still too high to see the sign changing of the oscil¬ 
lation pattern because the 7 band is minor where the 
expected sign changing temperature theoretically^ and 
experimentally^ usually located at 0.1T C ~ 150mK must 
be lowered. 

This kind of the band-dependent nodal structure dif¬ 
fers from the idea based on the symmetry protected nodal 
structure where the all plural bands are governed by the 
same gap symmetry. From this point of view, it is pos¬ 
sible that the d x 2_ y 2 nodal structure may not be truly 
sign-changing symmetry, rather an extended s-wave type 
with the anisotropic gap whose minima are along the 
( 110 ) direction. 

Some authors^ assert that the horizontal line nodes 
compatible to a chiral triplet state ( p x +ipy ) cos p z . From 
the present point of view, the Sommerfeld coefficient 
7 (H) behavior for H || c may not be consistent with 
this which gives a too larger 7 (H) at low field region 
than the experimental data— shown in Fig. 0] In this 
connection we should mention that no one succeeded in 
explaining the 7 (H) behavior— for H || ab. In particu¬ 
lar, it is viewed that the initial rise of H || ab is often 
assigned to the a + /? DOS because the plateau of 7 (H) 
for this mid field region H ~ 0.4 — 0.5T seems to corre¬ 
spond to a+ /3 DOS (43% of the total), which was taken 
as supporting evidence of the 7 major scenario. Since we 
cannot accept this view anymore, a full understudying of 
the 7 (H) behavior for H || ab remains mystery. 


G. Perspectives-Future experiments; Knight 

shift and FFLO 

There remain certainly important experiments to put 
forth the research front of this interesting material: NMR 
experiments have been extensively donoii, showing the 
absence of any change of the Knight shift— below T c . 
This is interpreted as freely rotatable d-vector in a spin 
triplet pairing, keeping always it perpendicular to an ex¬ 
ternal field which is as small as ~mT, a big mystery be¬ 
cause according to the recent ARPES^ the spin-orbit 
coupling which locks the d-vector to the crystalline lat¬ 
tice is an order of 200mV. Since the recent magnetization 
experiment— is directly able to detect the spin suscep¬ 
tibility change at the first order transition H c2a i, which 
amounts to ~ 10 % drop compared with the normal value. 
This is in sharp contrast with the Knight shift experiment 
by NMR which calls for reexamination of NMR experi¬ 
ments. 

Here we propose the Xj measurement to detect the 
FFLO stat o 49 : 50 expected to exist along H c2t ab below 
T ~0.8K where the first order transition is observed. Re¬ 
cently anomalous Tj _1 enhancement is observed^ when 
entering the FFLO state due to the appearance of the 
zero energy state at the domain walls where the FFLO 
order parameter is 7 r-shifted. There is a good chance to 
observe it if the second phase below H c2a b really exists, 
which we believe so. A necessary condition is that T 2 is 
short enough so that the spin-lattice relaxation Tj pro¬ 
cess is dominated through those zero energy states as 
in k-(BEDT-TTF) 2 Cu(NCS )2 case— In this connection 
we mention a recent /iSR experiment— that probes the 
peculiar vortex morphology at low fields of H || c and re¬ 
lated theory based on hidden criticality associated with 
multi-bandness— 


VII. SUMMARY AND CONCLUSION 

Based on the microscopic Eilenberger theory extended 
to a multiband case, we have studied the mixed state 
properties of a uniaxial anisotropic type II superconduc¬ 
tor, focusing on the interplay between the Pauli param¬ 
agnetic effect and multiband effect. A two band model 
calculation is set up and applied to Sr 2 Ru 04 - We have 
succeeded in reproducing the data both of the form fac¬ 
tors of SANS experiment s 12 : 27 and magnetic torque ex¬ 
periment— That leads us to the conclusion that to un¬ 
derstand the physics in S^RuCK it is indispensable to 
consider both the Pauli paramagnetic effect and multi¬ 
band effect simultaneously, which conspire to give rise to 
a variety of mysteries associated with the pairing sym¬ 
metry determination in this material. 

As agreed with the previous identification based on the 
single band analysis^, the pairing symmetry in S^RuCH 
is either singlet which is most likely or triplet with the 
d-vector locked in the ab plane which is less likely. The 
/3 ( 7 ) band is major (minor) with the mass anisotropy 60 
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(180). Namely, the /3 ( 7 ) band has a lager (smaller) gap. 
The gap structure is a full gap in the /3 band while in the 
7 band it is d x 2_ y 2 like. This simple picture was diffi¬ 
cult to reach because of the extreme two-dimensionality 
of this material which prevents conventional experimen¬ 
tal access. Now the dedicated and refined experimental 
tool o 12 i 28 i 52 i 53 which are able to align the magnetic field 
direction accurately within 1 ° enable us to uncover the 
physics of S^RuO-p 
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